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ABSTRACT 

We report the measurement of the spin-orbit angle of the extra-solar planets HAT-P-8 b, HAT-P-9 b, HAT-P-16 b and HAT-P-23 b, 
thanks to spectroscopic observations performed at the Observatoire de Haute-Provence with the SOPHIE spectrograph on the 1.93-m 
telescope. Radial velocity measurements of the Rossiter-McLaughlin effect show the detection of an apparent prograde, aligned orbit 
for all systems. The projected spin-orbit angles are found to be A =-lT'^1f^, -16° ± 8°, -10° ± 16°, -1-15° ± 22° for HAT-P-8, HAT-P-9, 
HAT-P-16 and HAT-P-23, respectively, with corresponding projected rotational velocities of 14.5 ± 0.8, 12.5 ± 1.8, 3.9 ± 0.8, and 7.8 
±1.6 kms"'. These new results increase to 37 the number of accurately measured spin-orbit angles in transiting extrasolar systems. 
We conclude by drawing a tentative picture of the global behaviour of orbital alignement, involving the complexity and diversity of 
possible mechanisms. 

Key words. Planetary systems - Techniques: radial velocities - Techniques: photometry - Stars: individual: HAT-P-8, HAT-P-9, 
HAT-P-16, HAT-P-23 



1. Introduction 

The number of discovered extrasolar planetary systems is con- 
stantly increasing, at a rate larger than 50 per year since 2007 
The sub-sample of confirmed transiting systems (~ 100) 
is explored with greater care, since the measurement of their 
real mass and radius allows studying their internal structure. 
For a sub-group of transiting systems, several teams have now 
collected the observation of spectroscopic transits with high- 
precision radial-velocity instruments, giving access to the rel- 
ative angle between the stellar spin plane and the planetary or- 
bital plane through the modeling of the Rossiter-McLaughlin ef- 
fect (Rossiter 1924). Unexpectedly, one third of the short-period 
planets for which those measurements are available have signif- 
icantly misaligned orbits with respect to their star's rotational 
axis, even polar or retrograde (e.g. |Triaud etal. (2010l; Hebrard 
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[erar('20T0l); |Winn et aL] ( |2010a| l). Such a distribution of stellar 
obliquities requires a scenario of formation and orbital evolution 
of extrasolar giant planets more diverse than pure one-planet mi- 
gration inside a proto-planetary disk. 

In this paper, we report the Rossiter-McLaughlin measure- 
ment of four additional transiting planets, HAT-P-8 b, HAT- 
P-9 b, HAT-P-16 b, HAT-P-23 b, obtained at Observatoire de 
Haute Provence with the SOPHIE high-precision radial-velocity 
(RV) spectrograph. The four systems discussed here are transit- 
ing extrasolar planets found by the HATNet photometric survey 
( [Bakos et al.|2002) 2004). 

HAT-P-8 b is an inflated hot Jupiter, with a mass Mp- 1.52 
Mjup, a radius Rp- 1.50 Rjup and a circular orbit of period 3.076 



days (Latham et al. 2009 j ). The host star HAT-P-8 is a solar- 
metallicity F main-sequence star with a significant projected ve- 
locity of about 12 km/s, and magnitude V =10.17. 

HAT-P-9 is an inflated, moderate-mass Jupiter-like planet, 
with Mp=0.78 Mjup and Rp- IA Rjup. It orbits a l ate F star in 
a 3.92289-day circular orbit ( [Shporer et "aL][2009| . The star is 
relatively faint, with V = 12.3. Ephemeris have been recently 
revisited (Dit tmann et al.|2010 i; we used the updated value for 
the observation scheduling and transit modeling. 
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HAT-P-16 b is a 4 Mjup planet transiting a V =10.8 mag- 
nitude star; its orbital period is 2.77596 days, its radius is 1.29 
Rjiip and it has a slightly eccentric orbit despite the short period 
(e=0.036 with a lOcr significance). The parent star has a mass of 
1.2 Mq and an effective temperature of 6158K (Buchhave et al 
[20T0I 1. 

HAT-P-23 b is an inflated and massive giant planet orbiting 
a GO dwaif star with V = 11.94 , in 1.212884d. Its mass and 
radius are 2.09 Mjup and 1.37 Rjup; although the orbital period is 
extremely short, the orbit of HAT-P-23 b is slightly but signifi- 
cantly eccentric, with e =0.106. Its detection and analysis have 
been reported by Bako s et al.| ( |20l0| l. 

Section|2]describes the spectroscopic observations, section 3 
presents additional photometry that was performed and analysed 
to refine the ephemeris of HAT-P-8 b, section |4] discusses the 
modeling of the Rossiter-McLaughlin effect. In section|5]we put 
these discoveries into context and discuss the global picture of 
orbital alignment. 



2. Spectroscopic observations 

The spectroscopic observations have been obtained at 
Observatoire de Haute Provence with the 1.93m telescope, in 
the framework of the large program (reference lOA.PNP.CONS) 



Table 1. Updated orbit and planetary parameters for HAT-P-8 b. 
A circular orbit is assumed. 



led by the SOPHIE Consortium, presented in detail by Bouchy 
[eTaTi'SOOga). The SOPHIE insti'ument (Perru chot et al.||2008| l 
is a fiber-fed environmentally stabilized echelle spectrograph 
covering the visible range from 387 to 694 nm. The spectral 
resolving power is 70,000 in the High-Resolution mode, and 
40,000 in the High-Efficiency mode. For high-cadence obser- 
vations required in a Rossiter sequence, the limit between both 
modes is around V - 9.5 - 10: the gain in resolution does not 
compensate the increase of readout noise for fainter stars. 

The observing template "objAB" is used, in order to simul- 
taneously monitor the sky background contamination. The cor- 
rection is required when the sky brightness is large and the Earth 
velocity is not sufficiently separated from the stellar velocity. 

The radial velocities are obtained by cross-correlating the 2D 
extracted spectra with a stellar numerical mask, following the 
method developed by|Baranne et al.|([l996[). 



2.1. HAT-P-8 

The transit of HAT-P-8 b has been observed during the night of 
19 July 2010. We observed the V = 10.17 mag HAT-P-8 star 
in the High-Efficiency mode, and slow readout mode to reduce 
the detector noise. Three measurements of this target were col- 
lected in 2009, and three in 2010 outside the transit night. These 
few measurements permitted to combine with a properly defined 
offset the orbital radial-velocity variations to the set of higher- 
accuracy Keck/HIRES data described in [Latham et al. | ( [2009| l. 
During the transit night (around JD=55397.5), we secured 35 
measurements, between 21:40 to 02:40 UT. Their signal-to- 
noise ratios range from 37 to 42, with exposure times from 220 
to 640 seconds. Observing at constant signal-to-noise ratio with 
SOPHIE allows reducing the radial-velocity systematics, that are 
due to the illumination-sensitive Charge Transmission Efficiency 
( [Bouchy et al. 2009b I. Although this systematic effect can be 
efficiently corrected after the proper calibration effect has been 
quantified, by observing at constant signal-to-noise ratio we re- 
duce the effect. The individual error bars range from 16 to 21 
m/s. The six bluest orders were removed for the velocity cal- 
culation, since they only add noise, and the G2 mask was used 



Parameters 


Values and l-cr error bars 


Unit 


Vr (SOPHIE) 


-22.383 ± 0.01 


kms"' 


P 


3.0763373 ± 0.0000031 


days 


K 


158.2 ±3.2 


ms"' 


Mp 


1.34±0.05 


Mj„p 


a 


0.0449 ± 0.0007 


AU 


T, (primary transit) 


2454437.67505 ± 0.00042 


BJD 



for the cross-correlation. The l-cr error bar of our measurements 
within transit is 18 m/s. It is similar for the measurements col- 
lected outside transit. 

The RV data of HAT-P-8 presented in this paper are shown 
on Figure [T] t ogether with the previous Keck/HIRES data of 
Latham et al.| (|2009|. They are phased to the ephemeris re- 
fined by our photometric observations (see section 3). Modeling 
the full out-of-transit radial-velocity sequence with both Keck 
and SOPHIE data with a Keplerian solution, we find a semi- 
amplitude of the variation 1.5-cr larger than in 'Latham et al. 
(2009), and values for the planetary mass and the semi-major 
axis that are varied accordingly. We update the orbital parame- 
ters and planetary mass in Table [T| The eccentricity being com- 
patible with at less than 2-cr level, we assume here a circular 



orbit, as in Latham et al. (2009). The bisector span is constant to 
the 2-0" level during the sequence. 

The velocity anomaly during the transit has an amplitude 
only 30% smaller than the orbital amplitude; this was expected 
from the significant projected velocity of the star and the large 
planetary radius (1 1.5 km/s and 1.5 Rjup respectively, in Latham 
|et al.[ ( |2009l l). This v sin /, can also be estimated from the cross- 
correlation function, as described in [Boisse et al.[ ( |20T0j l. We find 
V sin is- 12.9 + 1 km/s, larger but in agreement with the value re- 
ported in Latham et al. ( |2009) . 



2.2. HAT-P-9 

Two partial transits of HAT-P-9 have been recorded, on 28 
December 2010 and 1 January 201 1 in the High-Efficiency mode 
and slow readout. Ten (seventeen) measurements were secured 
during the first (resp. second) night. Exposure times ranging 
from 830 to 1200 seconds resulted in spectra with signal-to- 
noise ratio ranging from 14 to 31. An error of 12 m/s was 
quadratically added to the photon noise uncertainty, to account 
for instrumental systematics at low signal-to-noise ratio and get 
a final value close to 1 for the orbital fitting. We added the fif- 
teen measurements obtained with SOPHIE in 2007-2008, at the 
planet's discovery by [Shporer et al. ( [2009| l. All spectra were re- 
duced in the same way, using the FO correlation mask and spec- 
trograph orders 3 to 34. The values are reported in Table |2] We 
finally adjusted the planetary orbit using the last ephemeris de- 
rived by additional photometry by ,Dittmann et al. ( j2010| ). All 
RV data are shown on Figure [T] 

The transit-center epochs, using this ephemeris, are expected 
at Julian dates 5559.44531 and 245563.37109. Due to variations 
of the weather conditions, we had to interrupt the first sequence 
30 minutes before the end of the transit, and we did not get data 
after transit on the 28 December 2010 (observations lasted from 
20:17 to 23:54 UT). The second sequence nicely completes the 
first one, with a full transit observed and seven measurements af- 
ter transit (from 19:20 to 00:59 UT). No data were obtained be- 
fore transit however, because of the twilight. The radial-velocity 
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Fig. 1. SOPHIE (black plain circles), Keck/HIRES (red dia- 
monds), FIES (open green circles) radial-velocity data of HAT- 
P-8, HAT-P-9, HAT-P-1 6, and HAT-P-23. The superimposed or- 
bital solution is our best-fit to the data (not including the spec- 
troscopic transits in this plot). 



time series of both transits were assembled by letting a free off- 
set between both sequences, to allow for an observed shift of 37 
m/s. The reason for this may be related to the atmospheric chro- 
matic dispersor setup. This, however, does not cause additional 
jitter during the short timescale of the transit itself, thus all data 
are kept in the fit. The bisector span is constant at the 1.6-cr level 
during the sequence. The average width of the cross-correlation 
function is 16.6 km/s, which corresponds to a projected rota- 
tional velocity of 12.2 km/s, with a precision of about 1 km/s. 



This value is compatible with the one found by Shporer et al. 
(|2009jl from the analysis of stellar pai^ameters (1 1.9 ± 1 km/s) . 



|2010[ ). The bisector span is constant within 1 .4-cr during the se- 
quence. 

A projected stellar rotational velocity of 3.0 + 1 km/s is 
found from fitting the cross-correlation function, in agreement 
with the value reported in iBuchhave et alT (2010 1 (3.5 + 0.5 
km/s)). 

2.4. HAT-P-23 

The transit of HAT-P-23 b has been observed during the night 
of 25 August 2010, a few days after its announcement by Bakos] 
et aL] ( |201 1 1. In this sequence, SOPHIE was used in the High- 
Efficiency mode and fast readout mode. In order to maintain the 
signal-to-noise ratio constant (~ 24 ± 1) over the 6.3-hour se- 
quence, an exposure time of 787 to 1200 s has been used. A RV 
precision of 24 to 36 m/s has been obtained, after cross corre- 
lating with the G2 mask, including all orders. Due to the full 
Moon, there is a significant background contamination entering 
the spectrograph fibers. The stellar spectrum is a posteriori cor- 
rected for this effect, using the neighbour fiber monitoring the 
sky flux simultaneously to each exposure. The value of the cor- 
rection varies from 100 to 500 m/s during the sequence, increas- 
ing towards the end of the sequence when the moon light became 
more prominent. The velocity difference between the barycentric 
Earth RV and the target's RV is about 5 km/s during the tran- 
sit night. This contamination by the background light induces a 
residual jitter in the last data points of the series, but does not 
prevent the detection of the RM anomaly, as seen on Figure [T| 
The sequence on 25 August lasted from 20:01 to 02:20 UT, with 
a transit expected from 21:41 to 23:51 UT following ephemeris 
of 'Bakos et al. (2010 1. Twenty-three measurements were ob- 
tained during the transit night, plus one additional measurement 
two nights before the transit. The amplitude of the RM anomaly 
is of the order of lOOm/s and the orbit is prograde. Keck/HIRES 
radial velocities were added to our data set in order to simulta- 
neously fit the orbit and the transit (Figure [T] and section |4|i. The 
bisector span is constant within Icr during the sequence. 

We measure a projected stellar rotational velocity of 10 + 1 
km/s, larger by 2-cr than the value of 8. 1 ± 0.5 km/s reported in 
[Bakos et al.| ( pOTOl l. 



2.3. HAT-P-1 6 



3. Additional photometry 



The transit of HAT-P-1 6 b has been observed during the night 
of 15 August 2010 with SOPHIE in the High-Resolution mode 
and fast readout mode. In order to maintain the signal-to-noise 
ratio constant (~ 27 ± 1) over the 5-hour sequence, an expo- 
sure time of 848 to 1400 s has been used. The cross-correlation 
is calculated with a G2 mask which lowers the individual error 
bars, excluding the 14 bluest orders to diminish the systematic 
effects due to an incorrect positioning of the atmospheric chro- 
matic dispersor during this observation. A RV precision of 10 to 
15 m/s has been obtained. The stellar spectrum was not affected 
by any significant background light. One additional measure- 
ment was obtained the night before the transit in order to assess 
the RV uncertainty. Eighteen measurements were secured during 
the transit night. The sequence lasted from 22:19 to 03:27 UT. 
The transit was expected to occur between 23:47 and 02:51 UT, 
following ephemeris in |Buchhave et al. ( 2010[ l.The amplitude of 
the RV anomaly is extremely small, and the pattern appears reg- 
ular as for a prograde aligned orbit, with a redshift followed by 
a blueshift. Figure [T] shows the SOPHIE data overplotted with 
literature data obtained with FIES and HIRES ( |Buchhave et al.| 



The observation of the spectroscopic transit of HAT-P-8 b mo- 
tivated the scheduling of photometric observations of the tran- 
sit, since the egress seemingly occurred several minutes before 
expectations on the 19 July 2010 sequence according to the 
ephemeris by .Latham et al.| ( f2009 ). Photometric observations to 
refine the ephemeris were therefore carried out from several sites 
in France and Spain on 25, 28 and 31 August, from both ama- 
teur and professional observatories, and are briefly described in 
Table |5] Each lightcurve has been acquired and derived from 
observations independently; aperture photometry has been per- 
formed with optimized parameters for each system and condi- 
tions. Additionally, a transit light curve obtained on the night of 
September 13th by Ken Hose, was recovered from the Extrasolar 
Transit Database (ETD). 

The conversion from UTC timestamps to Barycentric 
Dynamical Time (TDB) standard was done using the code by 
Eastman et al.| ( [20T0| . Although not mentioned explicitly in the 
discovery paper, the KeplerCam timestamps were assumed to be 
provided in UTC and were therefore corrected to TDB by adding 
65.184 seconds to the measured central times (see |Eastman et al.] 
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|2010| l. Anyway, this correction does not change the ephemeris 
above the 1-cr level. 

Each light curve was fitted individually with the orbital and 
planetary parameters fixed to those by Latham et al. ( 2009| l in 
order to obtain the central times of transit. The error bars were 
estimated by means of the prayer bead method (see, for exam- 
ple. Desert et al.| ([2009 )), which consists in fitting a set of syn- 
thetic datasets obtained by shifting sequentially the residuals of 
the best-fit model rigidly and adding them back to the best-fit 
trapezoid function. The dispersion in the set of fit parameters 
obtained in this way is then used as an estimate of the uncer- 
tainties. This method permits to preserve the structure of the 
residuals and is therefore supposed to produce error estimates 
that take into account the covariant noise in the light curves. At 
each iteration of the processes, the limb-darkening coefficient 
corresponding to each photometric band in which a transit was 
observed was drawn randomly from a normal distribution cen- 
tred at the value obtained by interpolating the tables of Claret 
or 



(2000 



^taU 



Claret ( 2004 1 to the stellar parameters from Latham 
2009[ ). The dispersion of this distribution is defined by the 



extremes values of the limb-darkening coefficients obtained by 
allowing the stellar parameters to vary within the reported l-cr 
intervals. The KeplerCam lightcurves from the discovery paper 
( Latham et al.|2009 1 were fitted in a similar way. Figure|2]shows 
the light curves obtained by the different observers, together with 
the light curves from KeplerCam, binned to 10 minutes. 

An improved ephemeris was obtained by linear regression 
of the complete set of central times obtained, excluding the 
KeplerCam transit of November 1st 2007, since only the ingress 
was observed: 



To[BJD - 2450000] = 4437.67505 ± 0.00042 
P = 3.0763373 + 0.0000031 days 



(1) 
(2) 



with a covariance cov(f, To) - -5.1443 x 10"'*', and reduced 
= 0.75. The value of^^ probably means that the error bars 
of the central times of transits have been slightly overestimated. 
However, we decided not to scale them and therefore to report 
conservative ephemeris uncertainties. 

This ephemeris is in good agreement with those obtained by 
[Simpson et al.| ( [20T0l ), but the orbital period differs significantly 
(by almost 8-cr) from the value reported by |Latham et al. ( 2009| l. 
As a consequence, the predictions for the July 19th transit dif- 
fer by more than 19 minutes (8-cr). Using our computation, the 
expected central time of transit observed by SOPHIE is 

T,[BJD - 2450000] = 5397.49230 ± 0.00089 , 

where the covariance has been considered for propagating the 
error from the ephemeris. 

To explore the discrepancy with the prediction by [Latham 



et al. (2009 1, we adjusted the data of this target obtained by the 
HATnet telescopes together with the KeplerCam transits, in an 
attempt to reproduce the ephemeris reported in the discovery pa- 
per The orbital period thus obtained is in agreement with our 
determination to the 2-cr level and it differs from the period by 



Latham et al. ( 2009 1 by more than 7 cr. We conclude there exists 
a problem -probably a misprint- with the value reported in Table 
13 of [Latham et al. (2009 1, and that the observed discrepancy is 
not rea0 

^ Curiously, the orbital period reported on Table 13 of the version 
of the paper found on astro-ph lhttp://arxiv.org/pdf/0812.1 161 1 dilfers 
from the published value. The astro-ph version of the orbital period is 
closer to our determination (they differ by 3.4 cr) and the difference in 



Using the above values for To and P we adjusted the phase- 
folded curve using the analytical formulae of 'Mandel & Agol 
(2002) with a linear limb-darkening law to obtain the radius ra- 
tio k - Rp/Rt, the system scale a/R^ and the impact parame- 
ter b = a/Rf cos(i). The minimum value of the;^f^ statistic was 
obtained using the down-hill simplex method dNelder & Mead 
1965), and the errors were estimated using the Prayer Bead 
method. At each iteration of the processes, the limb-darkening 
coefficient corresponding to each photometric band in which a 
transit was observed was drawn randomly from a normal distri- 
bution centred at the value obtained by interpolating the tables 
oflCTaret (2000) or Claret (2004) to the stellar parameters from 
Latham et al.,(,2009) . The dispersion of this distribution is de- 
fined by the extreme values of the limb-darkening coefficients 
obtained by allowing the stellar parameters to vary within the 
reported 1-cr intervals. 

The obtained parameters are in perfect agreement with those 
from [Latham et al.| ( [2009 | l, although with a larger uncertainty. 
We thus use the transit parameters of the latter paper to model 
the spectroscopic transit. 



4. Modeling the spectroscopic transits 

Modeling of the radial-velocity anomaly has been performed us- 
ing the analytical approach developed by O hta et al. (200 5 ). The 
complete model has 10 parameters: the four standard orbital pa- 
rameters for circular orbits, the radius ratio Rp/Rt, the orbital 
semi-major axis to stellar radius a/Rt (constrained by the transit 
duration), the sky-projected angle between the stellar spin axis 
and the planetary orbital axis A, the sky-projected stellar rota- 
tional velocity v sin is, the orbital inclination /, and linear limb- 
darkening coefficients based on |Claret[ ( [200 4) tables for filter g' . 
Our goal here is to measure the spin-orbit projected angle A from 
the fits of the RM anomaly shape; we thus let free this parameter 
A in the fits. As in the case of HAT-P-6 ( Hebrard et al.|[20Tll l, 
the parameters having a significant impact on the RM fits and 
the A value are mainly the projected stellar velocity v sin is and 
the systemic velocity y of the data, and less critically, the RV 
semi-amplitude K and orbital inclination /. The K and / param- 
eters were fixed to their mean value in the fit but we explored 
the range of their error bars and quadratically added their small 
contribution to the final errors on v sin /^and A. We constructed a 
grid scanning the different values for A, v sin and the systemic 
velocity y, and computed the for each model. We performed 
a conservative re-scaling of the array, introducing a correc- 
tion for the x^ fluctuations due to the limited number of degrees 

~ 2002 ). Adopting the 

mmimum;^^^ 



of freedom, as explained in Hebrard et al 

obtained for each value on the grid of a given pa 



rameter (e.g. A), we determined the uncertainty on this parameter 
from the;\f^, Ax^ variations of 1, 4, and 9 determining the inter- 
val confidences at 1,2, and 3 cr. 

In RM fits, the degeneracy between vsin/, and A increases 
with the orbital inclination /. For impact parameters b that are 
not too small this correlation is taken into account in this work, 
by exploring the entire space of possible values for v sin /, and 
A. Figure [3[ shows a zoomed view of the spectroscopic transits 
and their best-fit solutions. The confidence interval contours esti- 
mated from rescaled-;^'^ variations against A and v sin / , are plot- 
ted inFig.[4[ 

the corresponding predictions of the July 19th 2010 is about 6.6 minutes 
(less than 3 cr). This reinforces our conclusion that the discrepancy is 
not significant, but stems rather from a typographical error. 
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Fig. 2. Upper panel: Transit light curves of HAT-P-8 in three dif- 
ferent bands: V (top; green), R (middle; red) and Sloan-z (bot- 
tom; brown), shifted veitically to improve legibility. Different 
symbols indicate different observers or nights, as indicated in 
the legend. The dashed curve is the best fit model for each pho- 
tometric band. Bottom panel: Residuals to the best fit model. The 
data has been shifted to avoid crowding. 



4.1. HAT-P-8 

We used the orbital parameters as derived in Section 2.1 and 
3 and the transit parameters from Latham et al. ( 2009| l. Two 
datasets were used in search for the best solution: first, the 
SOPHIE data set (35 measurements during the transit night), as 
described in Section 2.1; second, the data set composed of the 
SOPHIE data and the FIES measurements as recently published 
by jSimpson et al. ( 2010[ ). This latter study includes a total of 28 
data points. When quadratically adding 15 m/s to the FIES in- 
dividual errors (as in Simpson et al. ( |2010j l), and allowing for a 
free offset between both SOPHIE and FIES data sets, we obtain 
values of vsin/.,^ 14.5 + 0.8 km/s, A - -17° + 8 ° for the pro- 
jected rotational velocity and spin-orbit angle, respectively, with 
a reduced of 1.19. Fits with the extreme values of / show that 
an additional uncertainty should be quadratically added to the 
uncertainty obtained from the;^^^ variations. We checked that, as 
expected, variations of other parameters within their eiTor bars 
have no significant impact on A. The final value of A is - 




0,2 0,22 
H.4T-P- 16 



0,24 0,26 
phase 



0,2B 0-3 
SOPHIE 




The values reported by |Simpson et al. (20101 are -9.7° for A 
(l-o- range from -17.4° to -0.7°) and 1 1.8 + 0.5 km/s for v sin /j. 



phase 

Fig. 3. Close-in view of the radial-velocity measurements of 
HAT-P-8, HAT-P-9, HAT-P-16 and HAT-P-23 obtained with 
SOPHIE during transit nights. The best-fit Rossiter-McLaughlin 
models are superimposed as a red dashed line. FIES data of 
HAT-P-8 are included as green circles (Simpson et al. 2010[ l. 
Two separate observing nights are shown with different symbols 
for HAT-P-9 (filled circles: 28 December 2010, open circles: 1 
January 201 1). 



Adding the SOPHIE values therefore tends to slightly increase 
the spin-orbit angle (both estimates being in agreement within 
1-0"), while the derived vsin/, increases by about 2-cr. It is a 
known consequence in fitting spectroscopic transits of stars that 
are significantly rotationally broadened, that the v sin is is over- 
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estimated ( [Hirano et al.||2010) l, with consequences on the esti- 
mated A value that are of order of magnitude smaller than our 
error bar ( Triaud et al.|2009 ). Figure]?] shows the distribution of 
values against A for this system. 

4.2. HAT-P-9 

We first fit the planetary orbit, using the full set of SOPHIE data 
since 2007, and new ephemeris by Dittmann et al. ( 2010[ l: we 
find less dispersion in the residuals and a planetary mass slightly 
smaller than given in Shporer et al. (2009) . The RV best-fit semi- 
amplitude is K=73.8+ 6.7 m/s and Mp= 0.67 + 0.08 MjupWith 
M, = 1.28+0.13 Mq. 

For this planet, the semi-amplitude of the RM anomaly is 
larger than the amplitude of the orbital motion, a rare situation 
again explained by the large planetary radius, large stellar ro- 
tational velocity, and relatively light planet (hence small orbital 
amplitude). The Rossiter effect is then modeled by combining 
both December 2010 and January 20 11 sequences, and letting a 
free offset between these series. Half of the transit is observed 
twice, so that the offset is determined from a subset of about 
eight measurements per series. The offset is found to be 37 m/s 
in order to minimize the residuals. With the total sample of 27 
measurements during or close to the transit, the best-fit of the 
Rossiter-McLaughlin anomaly gives the following solution: ro 
tational velocity of 12.5 + 1.8 km/s (in agreement with jShporer 
eraL] ( ]2009| ) and spin-orbit angle /I of -16° + 8°. The reduced 



X^~^alue is 0.68. Taking into account the impact of RV semi- 
amplitude and inclination errors does not change significantly 
the value nor its error bar The transit data and its best fit are 
shown in Figure ]3] and the distribution of the derived quan- 



tity is shown on Figure |4] 



4.3. HAT-P-16 

Due to a lower value of the stellar rotational velocity and a non- 
inflated planet, the amplitude of the Rossiter anomaly on the star 
HAT-P-16 is small, while the individual measurement precision 
is greater compared to the other sequences. We perform the fit 
using the transit parameters as derived by Buchhave et al.| ( 2010| l 
and find a best-fit solution for a stellar projected rotational veloc- 
ity of 3.9 + 0.8 km/s (in agreement with Buchhave et al. ( |2010 1) 
and a spin-orbit angle of -10 ± 16 °. The reduced of :he fit is 
0.89, and the fitted data and;^'- distribution are shown on Figures 
]3]and]4] We checked that the variation of the RV semi-amplitude 
and inclination within their error bars did not have a noticeable 
impact on the value and error of the spin-orbit angle and rota- 
tional velocity in this case. Note that with an effective tempera- 
ture of 6158 K, one may expect the stellar rotational velocity to 
be greater than 3.9 km/s. An average value for such a star would 
rather be about 10 km/s ( Nordstrom et al.|2004| l. It is thus pos- 
sible that the stellar rotation axis of HAT-P-16 be inclined over 
the line of sight, resulting in a real spin-orbit angle larger than 
measured. 



4.4. HAT-P-23 

Finally, the SOPHIE data of HAT-P-23 were adjusted in the same 
way than previously, using the transit parameters and ephemeris 
of Bako s et al.| ( |2010| ). We derive a stellar rotational velocity 
of 7.8 ± 1.6 km/s and a projected angle between the orbital 
plane and the stellar equatorial plane, of 15° +22°. Again, we 
checked how the variations of the RV semi-amplitude and incli- 



nation within their error bars would affect the value and error of 
the spin-orbit angle and rotational velocity; here the error on the 
inclination has some impact, which is included in the final error 
bars. The vsin/^ derived from the fit is in good agreement with 
the value reported in Bakos et al. (2010 1 from stellar analyses. 
Figure ]4]shows the;^f distribution as a function of A. 

5. Discussion 

Table ]6] summarizes the values derived by our analyses of the 
spectroscopic transits of HAT-P-8 b, HAT-P-9 b, HAT-P-16 b, 
and HAT-P-23 b. The total number of systems with published 
projected spin-orbit angles is now 37. Ten of them are in the 
range \A \ - 20-120°, 7 are retrograde (more than ~120°), 22 are 
less than ~20°, among which the four systems presented here. 
Let us recall that the A value is the projected spin-orbit angle. 
The analysis of the real spin-orbit angle i// distribution by Triaud] 
et al. (20101 showed that 45 to 85 % of hot Jupiters may in fact 
have an orbital inclination greater than 30° . The recent study by 



Schlaufman ( 2010j l discusses the projected stellar rotation axis 
and the spin-orbit measurements and concludes that systems of 
massive stars and planets are more likely to show an apparent 
misalignment. 

Since the discovery of the first misaligned system ( [Hebrard 
|et al. 2008 ), there have been a lot of processes invoked for 



explaining them: three-body Kozai me chanism (Fabrycky & 
Tremaine 2007 1, planet-planet scattering ( Chatterjee et al.|2008 



Ford & Rasio||2008| l, tilted disks ( |Bate et al.||2010| |Lai et~aL 



20101, and insufficient tidal torque of stars hotter than 6200 K 



-hence a conservation of an original misalignement due to the 



preceeding mechanims for planets around hot stars ( | Winn et aT 



2010a I. Yet another process should be taken into account for 
this global picture: the tidal instability, which may affect all sys- 
tems by tilting the rotation axis of the star in relatively short 
timescales (Cebron et al. 201 1), and thus contribute to the distri- 
bution of orbital inclinations. The four systems analysed in this 
paper seem to be out of the so-called forbidden zone - which is 
defined by an orbital period less than 1/3 of the stellar rotational 
period (L e Bars et al. 2010) , and where the tidal instability can- 
not grow - and thus might have suffered a change of obliquity 
during their history. Complete simulations of the tidal evolution 
of the systems should be carried out in order to unveil possible 
impact of the tides onto the stars or the planets, particularly HAT- 
P-16 and HAT-P-23 which exhibit non-zero eccentricity despite 
their age (2.0 and 4.0 Gyr respectively). 

It is thus almost easier now to explain misaligned systems, 
than the more abundant hot Jupiters close to alignment. The ad- 
dition of four systems not strongly misaligned, as well as other 
recent discoveries renews our interest for this majority of sys- 
tems that are "surviving" the vast panel of possible causes for 
misalignment. Why should hot Jupiters get aligned with the 
equator of their star, and which parameter(s) would permit to 
predict the behaviour of the orbital planes? What constraint do 
we have on the sample of the aligned hot Jupiters, in terms of 
exclusion of other planets or stellar companions in the system, 
or stellar and disk properties? 

Figure ]6] (top) shows the distribution of the absolute value 
of the A projected spin-orbit angle with respect to the planetary 
mass, updated from a similar plot by |Hebrard et al. ( 2010j ) with 
almost a factor of two more systems. It shows that most of the 
systems show a low projected spin-orbit angle over the mass 
range of Jupiter-like planets. Among the planets in this mass 
range with misaligned orbits are included planets in binary sys- 
tems (as WASP-8, IQueloz et al.| ( |2010[ )) and seven systems in 
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retrograde orbits ( |/l | > 100°). These latter systems show some 
dependency with the planetary mass: only planets less massive 
than 3 Mjup are retrograde. 

There are at the moment: i) notable exceptions to any general 
rule (e.g. HAT-P-1 1 recently announced by | Winn et al. ( 2010b[ l, 
a A-103° orbit around a cool star, depicting an effect linked to 
the tidal dissipation rather than the sole effective temperature of 
the star) and ii) partial knowledge of the systems (unknown other 
planets or stellar companion in the system), that limit our un- 
derstanding of the common behaviour. Figure |6] (bottom) shows 
an updated version of the \A \ dependency with stellar effective 
temperature, as in Winn et al. (2010a). Our data might add a new 
exception of the rule that hot stars have planet with high obliq- 
uities (HAT-P-9 is has a low-angle misalignment, compatible at 
2-cr with zero with Te//^=6350K), although it may also indicate 
that the obliquity never changed during the system's evolution. 
Additionally, two of the systems reported in this paper lie close 
to the limiting temperature proposed by |Winn et al.| ( 2010b| l, 
HAT-P-8 and HAT-P-16 (with T,//=6200 and 6158 K respec- 
tively). Finally, HAT-P-23 with a cooler star (T,,//=5905K) and 
aligned orbit stands within the majority of such systems. The 
current sample of planets with measured spin-orbit projection 
angle hinders statistically significant conclusions concerning the 
difference between planets orbiting hot and cold host stars from 
being drawn. Indeed, considering T^eff - 6250 K as the limiting 
temperature between samples, a Kolmogorov-Smimov test on 
the data fails at rejecting the null hypothesis (i.e. that the sample 
of measured A for hot and cold systems are drawn from the same 
distribution) with a significance below a - O.ll^gQg, where the 
error bars are issued from 10000 bootstrap realisations of the 
data, using the uncertainties of effective temperature and A re- 
ported in the literature. This means that at the l-cr level, there 
exists a 28% of chance of erroneously rejecting the hypothesis 
that both populations are one and the same. Interestingly, we 
found slightly lower values for the significance when the limit- 
ing Teff is slightly higher: for Teff between 6300 K and 6450 K, 
the null hypothesis can be rejected at the l-cr level with a sig- 
nificance a - 0.2. We note, however, that for these temperatures 
only about 10 systems are classified as orbiting hot stars. 

The first steps are thus to systematically complete our un- 
derstanding of the systems, both by following up known tran- 
siting hot Jupiters in radial velocity on the long term, and by 
deriving detection limits of stellar companions by direct imag- 
ing. Companions able to efficiently misalign the planetary sys- 
tems through Kozai oscillations should generally be separated by 
no more than ~100 AU from the stars hosting the planets, oth- 
erwise the Kozai period exceeds that of the relativistic preces- 
sion, which then damps the Kozai cycles (Fabr ycky & Tremaine| 
2007| . Unfortunately, the transiting systems discovered during 
surveys such as HATnet or SuperWASP are commonly located 
at large distances (> 200 pc). Consequently, the adaptive optics 
imaging surveys, such as the ones carried out at Subaru ( [Narita 
et al.||2010) l or VLT (Ehrenreich et al., in prep.) have to probe 
the stellar environment within a few 0.1 arcsec. Searching for 
low-mass stars or brown dwarfs at these separations requires 
using innovative techniques such as angular differential imag- 
ing or sparse aperture masking (Lacour et al. 2011, in prep.). 
Alternatively, such companions could be detected through a 
long-term drift in the radial velocity of the star, providing they 
are observed long enough. A systematic collection of such data 
would provide useful constraints to the multi-body mechanisms 
and thus reinforce all other proposed scenarios for single planets 
around single stars. 
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Fig. 4. contour plot for the A values from the Rossiter- 
McLaughlin fits. 
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Table 5. Photometric observations of the transits of HAT-P-8 b used to refine the ephemeris. 



Observer 


Location 


Instrument 


Photometric Band 


Date of transit obs. 


Assoc. Adagio 


Toulouse, F 


82-cm, SBIG STL-6303E 


V 


08.28.10 and 08.31.10 


Hose 


Portland, Oregon 


32-cm, QSI516wsg 


R 


09.13.10 


Troncin 


OHP, F 


120-cm Tektronix TK1024 


V 


08.28.10 


Vanhuysse 


La Palma, S 


35-cm, SBIG STL-lOOle 


R 


08.28.10 



HAT-P-8b 




A = -17 



□ +9.2° 
-11.5' 



X = -lo-'iie- 




X = -16°±8° 


i 
i 


HAT-P-9b mk 


1 

i 

i 
i 




k- 

^ \ 

\ 
\ 

HAT-P-23b ^fll 


< 






\ 
\ 
\ 



Fig. 5. Sketch of the orbital and rotational orientations for all four systems discussed here, as viewed from Earth. 



Table 6. Spectroscopic transit modeled parameters of HAT-P-8, HAT-P-9, HAT-P-16, and HAT-P-23. Effective temperatures and 
planetary masses are taken from the discovery papers, except the planet mass of HAT-P-9 which is re-estimated in this work. 



Star 


VO 


V sin i , 


X 


2 

^red 




Mp 




km/s 


km/s 


degrees 




K 


Mj„p 


HAT-P-8 


-22.384 ± 0.0035 


14.5 ±0.8 


-17 (-11.5,-^9.2) 


1.19 


6200±80 


1.52±0.18 


HAT-P-9 


22.628 ± 0.008 


12.5 ±1.8 


-16± 8 


0.68 


6350±150 


0.67±0.08 


HAT-P-16 


-16.553 ± 0.0035 


3.9 ± 0.8 


-10.0 ± 16 


0.89 


6158±80 


4.193±0.094 


HAT-P-23 


-14.259 ± 0.007 


7.8 ± 1.6 


+ 15 ± 22 


1.02 


5905±80 


2.09±0.11 
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Fig. 6. Absolute value of the spin-orbit angle as a function of planetary mass (top) up to 25 Mjup and stellar effective temperature 
(bottom). The red filled large symbols show the new measurements. The uncertain values for CoRoT-lb, WASP-2 and Kepler-8 are 
included as green losenges. 
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Table 1. Radial velocity measurements obtained with SOPHIE 
ofHAT-P-8. 



JD-2,400,000. 


Radial Vel. 


Uncertainty 




[kms-i] 


[m s'M 


55072.45560 


-22.5463 


18.4 


55080.52188 


-22.2864 


18.5 


55088.49962 


-22.3445 


18.9 


55395.55227 


-22.4840 


20.8 


55396.52298 


-22.2275 


17.8 


55397.39997 


-22.3553 


18.5 


55397.40930 


-22.3904 


18.0 


55397.41822 


-22.3718 


18.6 


55397.43378 


-22.2944 


17.6 


55397.44022 


-22.2714 


17.5 


55397.44685 


-22.2634 


17.5 


55397.45327 


-22.2887 


18.3 


55397.46034 


-22.2989 


17.8 


55397.46719 


-22.2869 


17.2 


55397.47357 


-22.3169 


17.3 


55397.47896 


-22.2864 


17.1 


55397.48420 


-22.3135 


17.1 


55397.48938 


-22.3570 


17.5 


55397.49458 


-22.3382 


17.4 


55397.50028 


-22.3741 


17.3 


55397.51130 


-22.4156 


17.5 


55397.51675 


-22.4090 


17.5 


55397.52222 


-22.4740 


17.7 


55397.52780 


-22.4527 


17.8 


55397.53327 


-22.4539 


17.6 


55397.53856 


-22.4475 


17.8 


55397.54369 


-22.4959 


17.6 


55397.54868 


-22.4662 


17.7 


55397.55407 


-22.4575 


17.5 


55397.55914 


-22.4639 


17.6 


55397.56427 


-22.4359 


18.0 


55397.56930 


-22.4345 


18.1 


55397.57425 


-22.3988 


18.1 


55397.57918 


-22.3946 


18.4 


55397.58403 


-22.3932 


18.5 


55397.58884 


-22.3961 


18.3 


55397.59440 


-22.3973 


18.6 


55397.59935 


-22.3989 


18.5 


55397.60439 


-22.4194 


18.4 


55397.60947 


-22.4037 


18.4 


55398.53561 


-22.5182 


15.5 



Table 2. Radial velocity measurements obtained with SOPHIE 
ofHAT-P-9. 



JD-2,400,000. 


Radial Vel. 


Uncertainty 




[kms-'] 


[m s-i] 


55559.34192 


22.6812 


36.9 


55559.35866 


22.5790 


43.2 


55559.37506 


22.5852 


33.5 


55559.39135 


22.6719 


35.1 


55559.40763 


22.6588 


34.0 


55559.42391 


22.6997 


28.8 


55559.44366 


22.6173 


34.0 


55559.46000 


22.5507 


37.2 


55559.47627 


22.5503 


35.2 


55559.49255 


22.5119 


56.9 


55563.30211 


22.6347 


32.3 


55563.31818 


22.6929 


33.5 


55563.33464 


22.7358 


27.9 


55563.34868 


22.7057 


27.4 


55563.36333 


22.6758 


28.9 


55563.37770 


22.6792 


28.5 


55563.39233 


22.6012 


29.1 


55563.41194 


22.5295 


28.1 


55563.42434 


22.5684 


28.4 


55563.43747 


22.6301 


31.1 


55563.45125 


22.6167 


26.0 


55563.46384 


22.5940 


26.2 


55563.47781 


22.6040 


25.1 


55563.49242 


22.6422 


24.5 


55563.50598 


22.6353 


24.3 


55563.52041 


22.5892 


23.8 


55563.53780 


22.5987 


23.7 



Table 3. Radial velocity measurements obtained with SOPHIE 
ofHAT-P-16. 



JD-2,400,000. Radial Vel. Uncertainty 





[km s-i] 


[m s-'] 


55424.43002 


-16.4335 


10.7 


55424.44167 


-16.4491 


10.4 


55424.45587 


-16.4358 


12.5 


55424.47137 


-16.4391 


15.4 


55424.48567 


-16.5030 


15.4 


55424.50068 


-16.4913 


14.5 


55424.51498 


-16.5081 


11.6 


55424.52851 


-16.5096 


10.8 


55424.54198 


-16.5225 


11.2 


55424.55597 


-16.5680 


11.3 


55424.56858 


-16.5967 


11.0 


55424.57975 


-16.6180 


10.9 


55424.59027 


-16.6271 


10.9 


55424.60097 


-16.6502 


10.9 


55424.61203 


-16.6542 


11.1 


55424.62317 


-16.6594 


10.8 


55424.63377 


-16.6644 


10.9 


55424.64422 


-16.6681 


10.8 
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Table 4. Radial velocity measurements obtained with SOPHIE 
ofHAT-P-23. 



JD-2,400,000. 


Radial Vel. 


Uncertainty 




[kms-i] 


[m s'M 


55434.33411 


-14.0555 


24.5 


55434.34457 


-14.0666 


25.6 


55434.35558 


-14.0820 


24.5 


55434.36645 


-14.0898 


24.9 


55434.37731 


-14.1137 


25.2 


55434.38856 


-14.1572 


26.2 


55434.40002 


-14.2141 


26.1 


55434.41127 


-14.1490 


25.6 


55434.42196 


-14.1922 


25.3 


55434.43227 


-14.1831 


25.4 


55434.44326 


-14.2378 


25.8 


55434.45460 


-14.3113 


26.1 


55434.46798 


-14.4068 


27.1 


55434.47968 


-14.4386 


25.5 


55434.49150 


-14.4160 


25.5 


55434.50469 


-14.3966 


26.2 


55434.51904 


-14.4791 


27.4 


55434.53355 


-14.4316 


26.9 


55434.54732 


-14.5136 


25.9 


55434.55997 


-14.5404 


27.2 


55434.57210 


-14.4825 


24.5 


55434.58393 


-14.5221 


29.0 


55434.59768 


-14.5325 


35.7 



